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We present a family of methods, analytical and numerical, which can describe be- 
haviour in (non) equilibrium ensembles, both classical and quantum, especially in the 
complex systems, where the standard approaches cannot be applied. We demonstrate 
the creation of nontrivial (meta) stable states (patterns), localized, chaotic, entangled 
or decoherent, from basic localized modes in various collective models arising from the 
quantum hierarchy of Wigner-von Neumann-Moyal-Lindblad equations, which are the 
result of "wignerization" procedure of classical BBGKY hierarchy. We present the ex- 
plicit description of internal quantum dynamics by means of exact analytical / numerical 
computations. 

Keywords: Localization; pattern formation; multiscalcs; multiresolution; waveletons; 
(non) equilibrium ensembles. 



1. Localized Modes ("continuous qudits"): Why Need We Them? 

It is widely known that the currently available experimental techniques (and, appar- 
ently, those which will become avaiable in the nearest future) in the area of quantum 
physics as a whole and in that of quantum computations in particular, as well as 
the present level of understanding of phenomenological models, outstripped the ac- 
tual level of mathematical/theoretical description"'^. Considering, for example, the 
problem of describing the realizable states (Refs. [2H7]), one should not expect that 
planar waves and (squeezed) gaussian coherent states would be enough to character- 
ize such complex systems as qCPU (quantum Central Processor Unit)-like devices. 
Complexity of the set of relevant states, including entangled (chaotic) states is still 
far from being clearly understood and moreover from being realizable. As a starting 
point for our approach let us consider the following well-known example of GKP 
(Gottesman, Kitaev, Preskill)P scheme with DV (Discrete Variables) /qubit (with 
finite-dimensional code space embedded in the infinite-dimensional Hilbert space) 
or CV (Continuous Variables) for (optical) quantum computations, containing as 
a part (optical) nonlinearities, described by Kerr interaction or more general poly- 
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nomial Hamiltonians which are needed to reahze the state preparation and provide 
the process of CV quantum computatiorPI. It is an important example because: 

(a) its classical counterpart is described by polynomial Hamiltonians; 

(b) the proper qudits or building states (DV or CV) are well localized (but not 
well-defined mathematically, as we shall explain later). 

One of the questions which motivated our approach is whether it is possible to 
keep (a) and at the same time improve (b). Our other motivations arise from the 
following general questions: 

(A) How can we represent well localized and reasonable state in mathematically 
correct form? 

(B) Is it possible to create entangled and other relevant states by means of these 
new building blocks? 

In GKP scheme unphysical and not clearly defined mathematically logical qubit 
states are represented via infinite series of S functions: |0 >= X]^-oo '^(■^ ~ 
2sy^)\x >,|1 >= J27L-oo^{^ ~ (2s + > and approximated by the set 

of gaussian envelopes: < a:|0 >= Nq Y^tZ e-i/2(2'*'=v^)' e"^/^^ , < x\l >= 

Due to numerous mathematical and computational reasons, some of which are 
described below, such and related choices cannot be appropriate neither as a start- 
ing point on the route to the real qCPU device nor as a satisfactory theoretical 
description. So, it would appear that a first step in this direction is to find a reason- 
able extension of understanding of the quantum dynamics as a whole. One needs 
to sketch up the underlying ingredients of the theory (spaces of states, observables, 
measures, classes of smoothness, quantization set-up etc) in an attempt to provide 
the maximally extendable but at the same time really calculable and realizable de- 
scription of the dynamics of quantum world. The general idea is rather simple: it 
is well known that the idea of "symmetry" is the key ingredient of any reasonable 
physical theory from classical (in)finite dimensional (integrable) Hamiltonian dy- 
namics to different sub-planckian models based on strings (branes, orbifolds etc.) 
During the last century kinematical, dynamical and hidden symmetries played the 
key role in our understanding of physical process. Roughly speaking, the represen- 
tation theory of underlying symmetry (classical or quantum, groups or (bi)algebras, 
finite or infinite dimensional, continuous or discrete) is a proper instrument for de- 
scription of proper (orbital) dynamics. A starting point for us is a possible model 
for (continuous) "qudit" with subsequent description of the whole zoo of possible 
realizable (controllable) states/patterns which may be useful from the point of view 
of quantum experimentalists and engineers. The proper representation theory is 
well known as "local nonlinear harmonic analysis" , in particular case of simple un- 
derlying symmetry-affine group-aka wavelet analysis. From our point of view the 
advantages of such approach are as follows: 

i) natural realization of localized states in any proper functional realization of 
(Hilbert) space of states. 
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ii) hidden symmetry of chosen reahzation of proper functional model provides the 
(whole) spectrum of possible states via the so-called multiresolution decomposition. 

So, indeed, the hidden symmetry (non-abelian affine gronp in the simplest case) 
of the space of states via proper representation theory generates the physical spec- 
trum and this procedure depends on the choice of the functional realization of the 
space of states. It explicitly demonstrates that the structure and properties of the 
functional realization of the space of states are the natural properties of physical 
world at the same level of importance as a particular choice of Hamiltonian, or the 
equation of motion, or the action principle (variational method). At the next step 
we need to consider the consequences of our choice i), ii) for the algebra of observ- 
ables. In this direction one needs to mention the class of operators we are interested 
in to present proper description for a class of maximally generalized but reason- 
able class of problems. It seems that these must be pseudodifferential operators, 
especially if we underline that in the spirit of points i), ii) above we need to take 
Wigner-Weyl framework for constructing basic quantum equations of motions. It is 
obvious, that consideration of symbols of operators instead of operators themselves 
is the starting point as for the mathematical theory of pseudodifferential operators 
as for quantum dynamics formulated in the language of Wigner-like equations. It 
should be noted that in such picture we can naturally include the effects of self- 
interaction ( "quantum non-linearity" ) on the way of construction and subsequent 
analysis of nonlinear quantum models. So, our consideration will be in the frame- 
work of (Nonlinear) Pseudodifferential Dynamics {'^DOD). As a result of i), ii), 
we'll have: 

iii) most sparse, almost diagonal, representation for a wide class of operators in- 
cluded in the set-up of the whole problems. 

It's possible by using the so-called Fast Wavelet Transform representation for 
algebra of observables. 

Then points i)-iii) provide us by 

iv) natural (non-perturbative) multiscale decomposition for all dynamical quanti- 
ties, as states as observables. 

The simplest case we will have, obviously, in Wigner-Weyl representation. Ex- 
istence of such internal multiscales with different dynamics at each scale and tran- 
sitions, interactions, and intermittency between scales demonstrates that quantum 
mechanics, despite its linear structure, is really a serious part of physics from the 
mathematical point of view. It seems, that well-known underlying quantum com- 
plexity is a result of transition by means of (still rather unclear) procedure of quan- 
tization from complexity related to nonlinearity of classical counterpart to the rich 
pseudodifferential (more exactly, microlocal) structure on the quantum side. 

We divide all possible configurations related to possible solutions of our quantum 
equation of motion (Wigner-like equations, mostly) into two classes: 

(a) standard solutions; (b) controllable solutions (solutions with prescribed 
qualitative type of behaviour). 
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Anyway, the whole zoo of solutions consists of possible patterns, including very 
important ones from the point of view of underlying physics: 

v) localized modes (basis modes, eigenmodes) and constructed from them chaotic 
or entangled, decoherent (if we change Wigner equation for (master) Lindblad one) 
patterns. 

It should be noted that these bases modes are nonlinear in contrast with usual 
ones because they come from (non) abelian generic group while the usual Fourier 
(commutative) analysis starts from C/(l) abelian modes (plane waves). They are 
really "eigenmodes" but in sense of decomposition of representation of the under- 
lying hidden symmetry group which generates the multiresolution decomposition. 
The set of patterns is built from these modes by means of variational procedures 
more or less standard in mathematical physics. It allows to control the convergence 
from one side but, what is more important, 

vi) to consider the problem of the control of patterns (types of behaviour) on the 
level of reduced (variational) algebraical equations. 

We need to mention that it is possible to change the simplest generic group of 
hidden internal symmetry from the affine (translations and dilations) to much more 
general, but, in any case, this generic symmetry will produce the proper natural high 
localized eigenmodes, as well as the decomposition of the functional realization of 
space of states into the proper orbits; and all that allows to compute dynamical 
consequence of this procedure, i.e. pattern formation, and, as a result, to classify 
the whole spectrum of proper states. 

For practical reasons controllable patterns (with prescribed behaviour) are the 
most useful. We mention the so-called waveleton-like pattern which we regard as 
the most important one. We use the following allusion in the space of words: 

{waveleton}:={soliton} |J {wavelet} 

It means: 

vii) waveleton « (meta)stable localized (controllable) pattern 
To summarize, the approach described below allows one 

viii) to solve wide classes of general '^DOD problems, including generic for quantum 
physics Wigner-like equations, and 

ix) to present the analytical/numerical realization for physically interesting pat- 
terns. 

We would like to emphasize the effectiveness of numerical realization of this 
program (minimal complexity of calculations) as additional advantage. So, items 
i)-ix) point out all main features of our approach, Refs. [2H71 

2. Motivations 

2.1. Class of Models 

Here we describe a class of problems which can be analysed by methods described 
in Introduction. We start from individual dynamics and finish by (non)-equilibrium 
ensembles. All models belong to the "^DOD class and can be described by finite or 
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infinite (named hierarchies in such cases) system of ^DOD equations: 

a) . Individual classical/quantum mechanics (cM/qM): Hnear/nonlinear; 

{cM} C {qM}, * - quantized for the class of polynomial Hamiltonians 
H{p,q,t) = Y^i.j a,j{t)p'qi. 

b) . QFT-like models in framework of the second quantization (dynamics in 

Fock spaces). 

c. ) Classical (non) equilibrium ensembles via BBGKY Hierarchy (with reduc- 

tions to different forms of Vlasov-Maxwell/Poisson equations). 

d. ) Wignerization of a): Wigner-Moyal-Weyl-von Neumann-Lindblad. 

e. ) Wignerization of c): Quantum (Non) Equilibrium Ensembles. 

Important remarks: points a)-e) are considered in ^DO picture of (Non)Linear 
^DO Dynamics (surely, all qM C ^DOD); dynamical variables/observables are the 
symbols of operators or functions; in case of ensembles, the main set of dynamical 
variables consists of partitions (n-particle partition functions). 

2.2. Effects we are interested in 

(i) Hierarchy of internal/hidden scales (time, space, phase space). 

(ii) Non-perturbative multiscales: from slow to fast contributions, from the coarser 
to the finer level of resolution/decomposition. 

(iii) Coexistence of hierarchy of multiscale dynamics with transitions between scales. 

(iv) Realization of the key features of the complex quantum world such as the 
existence of chaotic and/or entangled states with possible destruction in 
"open/dissipative" regimes due to interactions with quantum/classical environ- 
ment and transition to decoherent states. 

At this level we may interpret the eflfect of mysterious entanglement or "quan- 
tum interaction" as a result of simple interscale interaction or intermittency (with 
allusion to hydrodynamics), i.e. the mixing of orbits generated by multiresolution 
representation of hidden underlying symmetry. Surely, the concrete realization of 
such a symmetry is a natural physical property of the physical model as well as 
the space of representation and its proper functional realization. So, instantaneous 
interactions (or transmission of "quantum bits" or "teleportation" ) materialize not 
in the physical space-time variety but in the space of representation of hidden 
symmetry along the orbits/scales constructed by proper representations. Dynam- 
ical/kinematical principles of usual space-time varieties, definitely, do not cover 
kinematics of internal quantum space of state or, in more weak formulation, we still 
have not such explicit relations. 

One additional important comment: as usual in modern physics, we have the 
hierarchy of underlying symmetries; so our internal symmetry of functional realiza- 
tion of space of states is really not more than kinematical, because much more rich 
algebraic structure, related to operator Cimtz algebra and quantum groups, is hid- 
den inside. The proper representations can generate much more interesting effects 
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than ones described above. We will consider it elsewhere but mention here only how 
it can be realized by the existing functorial maps between proper categories: 

{QMF} — > Loop groups — > Cuntz operator algebra — > Quantum Group 
structure, where {QMF} are the so-called quadratic mirror filters generating the 
realization of multiresolution decomposition/representation in any functional space; 
loop group is well known in many areas of physics, e.g. soliton theory, strings etc, 
roughly speaking, its algebra coincides with Virasoro algebra; Cuntz operator al- 
gebra is universal C* algebra generated by N elements with two relations between 
them; Quantum group structure (bialgebra, Hopf algebra, etc) is well known in 
many areas because of its universality. It should be noted the appearance of natural 
Fock structure inside this functorial sequence above with the creation operator re- 
alized as some generalization of Cuntz-Toeplitz isometries. Surely, all that can open 
a new vision of old problems and bring new possibilities. 

We finish this part by the following qualitative definitions of key objects (pat- 
terns). Their description and understanding in different physical models is our main 
goal in this direction. 

• By localized states (localized modes) we mean the building blocks for solutions 
or generating modes which are localized in maximally small region of the phase 
(as in c- as in q-case) space. 

• By an entangled/chaotic pattern we mean some solution (or asymptotics of solu- 
tion) which has random-like distributed energy (or information) spectrum in a full 
domain of definition. In quantum case we need to consider additional cntanglcd- 
like patterns, roughly speaking, which cannot be separated into pieces of sub- 
systems. 

• By a localized pattern (wavclcton) we mean (asymptotically) (meta) stable solu- 
tion localized in a relatively small region of the whole phase space (or a domain 
of definition) . In this case the energy is distributed during some time (sufficiently 
large) between only a few localized modes (from point 1). We believe it to be a 
good model for plasma in a fusion state (energy confinement) or a model for quan- 
tum continuous "qubit" or a result of the decoherence process in open quantum 
system when the full entangled state degenerates into localized (quasiclassical) 
pattern. 

2.3. Methods 

(i) Representation theory of internal/hidden/underlying symmetry, Kinematical, 
Dynamical, Hidden. 

(ii) Arena (space of representation): proper functional realization of (Hilbert) space 

of states. 

(iii) Harmonic analysis on (non)abelian group of internal symmetry. Lo- 
cal/Nonlinear (non-abelian) Harmonic Analysis (e.g, wavelet/gabor etc. analy- 
sis) instead of linear non-localized ?7(1) Fourier analysis. Multiresolution (mul- 
tiscale) representation. Dynamics on proper orbit/scale (inside the whole hi- 
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erarchy of multiscales) in functional space. The key ingredients are the ap- 
pearance of multiscales (orbits) and the existence of high-localized natural 
(eigen) modes'^. 

(iv) Variational formulation (control of convergence, reductions to algebraic sys- 
tems, control of type of behaviour). 

3. Set-up/Formulation 

Let us consider the following generic 5*000 dynamical problem 

L^iOp'j'f = 0, 

described by a finite or infinite number of equations which include general classes 
of operators Op* such as differential, integral, pseudodifferential etc 

Surely, all Wigner-like equations/hierarchies are inside. 

The main objects are: 

(i) (Hilbert) space of states, H = {5'}, with a proper functional realization, e.g.,: 
L2, Sobolev, Schwartz, C°, C^ ... C°°, Definitely, L^{R^), L^{S^), L^{S^ x 
5^), L'^{S^ X 5^ k Zn) are different objects proper for different physics inside. 

(ii) Class of smoothness. The proper choice determines natural consideration of 
dynamics with/ without Chaos/Fractality property. 

(iii) Decompositions 

i 

via high-localized bases (wavelet families, generic wavelet packets etc), frames, 
atomic decomposition (Fig. 1) with the following main properties: (exp) control 
of convergence, maximal rate of convergence for any 5* in any H. 

(iv) Observables/Operators (ODO, PDO, *D0, SIO,..., Microlocal analysis of 
Kashiwara-Shapira (with change from functions to sheafs)) satisfy the main 
property - the matrix representation in localized bases 

< ^'lOp'l^' > 

is maximum sparse: 

...\ 

D22 .. 

^?33 • • 

V : : ; 

This almost diagonal structure is provided by the so-called Fast Wavelet Trans- 
form. 

(v) Measures: multifractal wavelet measures {fii} together with the class of 
smoothness are very important for analysis of complicated analytical be- 
haviour. 
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Fig. 1. Localized modes. 

(vi) Variational/Projection methods, from Galerkin to Rabinowitz minimax, Floer 
(in symplectic case of Arnold- Weinstein curves with preservation of Pois- 
son/symplectic structures). Main advantages are the reduction to algebraic 
systems, which provides a tool for the smart subsequent control of behaviour 
and control of convergence. 

(vii) Multiresolution or niultiscalc decomposition, MRA (or wavelet microscope) 
consists of the understanding and choosing of 

1) . (internal) symmetry structure, e.g., affine group = {translations, dila- 

tions} or many others; construction of 

2) . representation/action of this symmetry on H = {5'}. 

As a result of such hidden coherence together with using point vi) we'll 
have: 

a). LOCALIZED BASES b). EXACT MULTISCALE DECOMPO- 
SITION with the best convergence properties and real evaluation of 
the rate of convergence via proper "multi-norms" . 

Figures 2, 3, 5, 6 demonstrate MRA decompositions for one- and multi-kicks 
while Figures 4 and 7 present the same for the case of the generic simple fractal 
model, Riemann-Weierstrass function [9]. 

(viii) Effectiveness of proper numerics: CPU-time, HDD-space, minimal complexity 
of algorithms, and (Shannon) entropy of calculations are provided by points 
i)-vii) above. 

(ix) Quantization via * star product or Deformation Quantization. 
The corresponding class of individual Hamiltonians has the form 

H{p,q) = ^ + Uip,q), (1) 
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where C/(p, q) is an arbitrary polynomial function on p, q, and plays the key role in 
many areas of physics Our starting point is the general point of view of a defor- 
mation quantization approach at least on the naive Moyal/Weyl/Wigner level. The 
main point of such approach is based on ideas from^l], which allow to consider the al- 
gebras of quantum observables as the deformations of commutative algebras of clas- 
sical observables (functions). So, if we have the classical counterpart of Hamiltonian 
(1) as a model for classical dynamics and the Poisson manifold Af (or symplectic 
manifold or Lie coalgebra, etc) as the corresponding phase space, then for quan- 
tum calculations we need first of all to find an associative (but non-commutative) 
star product * on the space of formal power series in h with coefficients in the 
space of smooth functions on M such that f * g = fg + fi{f, g} + X]n>2 5)i 
where {/, 5} is the Poisson brackets, Bn are bidifferential operators. In this paper 
we consider the calculations of the Wigner functions W{p, q, t) (WF) corresponding 
to the classical polynomial Hamiltonian H{p, q, t) as the solution of the Wigner-von 
Neumann equation 

ih^W *W ~W *H (2) 
dt ^ ' 

and related Wigner-like equations for different ensembles. According to the Weyl 
transform, a quantum state (wave function or density operator p) corresponds to the 
Wigner function, which is the analogue in some sense of classical phase-space distri- 
bution ^ . Wigner equation (2) is a result of the Weyl transform or "wignerization" 
of von Neumann equation for density matrix. 

Finally, such Variational-Multiscale approach based on points i)-ix) provides us 
by the full ZOO of PATTERNS: LOCALIZED, CHAOTIC/ENTANGLED, etc. 

In next Sections we will consider details for important cases of Wigner-like equa- 
tions. 

We present the explicit analytical construction for solutions of c- and q- 
hierarchies and their important reductions starting from quantization of c-BBGKY 
(Born-Bogolyubov-Green-Kirkwood-Yvon) hierarchy, which is based on tensor al- 
gebra extensions of multiresolution representation for states and observables and 
variational formulation. We give explicit representation for hierarchy of n-particle 
reduced distribution functions in the base of high-localized generalized coherent (re- 
garding underlying generic symmetry (affine group in the simplest case)) states given 
by polynomial tensor algebra of our basis functions (wavelet families, wavelet pack- 
ets), which takes into account contributions from all underlying hidden multiscales 
from the coarsest scale of resolution to the finest one to provide full information 
about (quantum) dynamical process. The difference between classical and quantum 
case is concentrated in the structure of the set of operators included in the set-up 
and, surely, depends on the method of quantization. But, in the naive Wigner- Weyl 
approach for quantum case the symbols of operators play the same role as usual 
functions in classical case. In some sense, our approach for ensembles (hierarchies) 
resembles Bogolyubov's one and related approaches but we don't use any perturba- 
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tion technique (like virial expansion) or linearization procedures. Most important, 
that numerical modeling in all cases shows the creation of different internal (coher- 
ent) structures from localized modes, which are related to stable (equilibrium) or 
unstable type of behaviour and corresponding pattern (waveletons) formation. 



4. BBGKY/Wigner Ensembles: from c- to q-cases 

We start from set-up for kinetic BBGKY hierarchy as c-couiiterpart of proper q- 
hierarchy. Let M be the phase space of ensemble of N particles (dimM = 6N) with 
coordinates x, = {qi.Pi), i = l,...,N,qi = {qj,qf,qf) G R^,p, = {p},pf,p^i) G 
R^^q = (^ij-'-j^w) G B?'^ . Individual and collective measures are: jii = Axi = 
dgijdpi, ^ = Hi^i Mi while distribution function Dn{xi, . . . ,XN;t) satisfies Li- 
ouville equation of motion for ensemble with Hamiltonian Hn and normalization 
constraint. Our constructions can be applied to the following general Hamiltonians: 

AT 2 

^^ = E(& + ^^(9))+ E U,j{q,,qj), (3) 
i=l i<i<j<N 

where potentials Ui{q) = Ui{qi, . . . , gjv) and Uij{qi,qj) are not more than rational 
functions on coordinates. Let Lg and Lij be the standard Liouvillean operators and 

Fs{xi,...,Xs;t) = V' I DN{xi,...,XN\t) /Zj (4) 

^ s+l<i<N 

be the hierarchy of reduced distribution functions, then, after standard manipula- 
tions, we arrive to c-BBGKY hierarchy: 

— + L^Fs = - I d/Xg+i ^ Li^s+iFs+\. (5) 
^ '"J i=i 

So, the proper dynamical formulation is reduced to the (infinite) set of equations for 
correlators/partition functions. Then by using physical motivated reductions or /and 
during the corresponding cut-off procedure we obtain, instead of linear and pseudod- 
ifferential (in general case) equations, their finite-dimensional but nonlinear approx- 
imations with the polynomial type of nonlinearities (more exactly, multilinearities) . 
To move from c- to g-case, let us start from the second quantized representation 
for an algebra of observables A = (Aq, ^i, . . . , Ag, ...) in the standard form A = 
Aq + ^ dxi^+{xi)Ai^{x-i) + . . . +(s!)-^ / dxi . . .dxs^+{xi) . . .^+{xs)A,'^{x,) . . . 
^{xi) + . . . N-particle Wigner functions 

Ws{xi,...,Xs) = I dki . . .dksex.p{- i'^kiPi)TTp'^+{qi - ^hki) . . . (6) 

•' 1=1 

- ^/ifcs)*(gs + ^hks) . . . *(gi ^hks) 
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allow us to consider them as some quasiprobabilities and provide useful bridge 
between c- and q-cases: 



oo . s 

< A >=TrpA = ^(s!)~i / '[[d^iiAs{xl,...,Xs)Ws{xl,...,Xs). 



(7) 



The full description for quantum ensemble can be done by the whole hierarchy of 
functions (symbols): W = {Ws{xi, . . . , Xg), s — 0, 1, 2 . . .} So, we may consider the 
following q-hierarchy as the result of "wignerization" procedure for c-BBGKY one: 

s s 
dtWs{t,Xl, ...,Xs) ^^L°Ws{xi, ...,Xs) + Y,Yl Lj,nWs{xi, ...,Xs) (8) 
j — 1 j<n n—1 



/ dxs+iS{ks+i)Lj^s+iWs+i{xi, . . . ,Xs+i), 



= -{^)kjpj, L,.n = {ih)-^ I MVi 



(9) 



In quantum statistics the ensemble properties are described by the density operator 
p{t)^Y.w^\^^{t)><-^,{t)\, ^u;, = 1. (10) 

i i 

After Weyl transform wc have the following decomposition via partial Wigner func- 
tions Wi{p,q,t) for the whole ensemble Wigner function: 



i 

where the partial Wigner functions 



(11) 



are solutions of proper Wigner equations: 



dt 



m dq 



E 



(2£+l)! 



dp 



(13) 



Our approach, presented below, in some sense has allusion on the analysis of 
the following standard simple model considered in Ref. [TJ Let us consider model 
of interaction of nonresonant atom with quantized electromagnetic field: H = 
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^ + U{x), U{x) = Uo{z,t)g{x)a^a, where potential U depends on cre- 
ation/annihilation operators and some polynomial on x operator function (or ap- 
proximation) g{x). It is possible to solve Schroedinger equation ihd\'i'>/dt ~ 
Hl'i' > by the simple ansatz 

\^{t) / dx'^„{x,t)\x > (8)|n >, (14) 

-oo 

which leads to the hierarchy of analogous equations with potentials created by n- 
particle Pock subspaces 

where 5'„(a;,i) is the probability amplitude of finding the atom at the time t at 
the position x and the field in the n Fock state. Instead of this, we may apply the 
Wigner approach starting with proper full density matrix |5'(t) >< ^(i)|: 

P=Y1 Wn'w:^,, J Ax' j Ax"-^n'{x' ,t)-^ln{:x" ,t)\x' ><x"\®\n' ><n"\. (16) 

n' ,n" 

Standard reduction gives pure atomic density matrix 



Pa EE / < n\p\n >= (17) 

J 71=0 

J2\^n\^ J J dx"^„ix',t)Kix",t)\x' ><X"\. 

Then we have incoherent superposition 

oo 

W{x,p, 

n=0 

of the atomic Wigner functions (12) corresponding to the atom motion in the po- 
tential Un{x) (which is not more than polynomial in x) generated by n-level Fock 
state. They are solutions of proper Wigner equations (13). The next case describes 
the important decoherence process. Let us have collective and environment subsys- 
tems with their own Hilbert spaces TC = Tic ® Ti-e Relevant dynamics is described 
by three parts including interaction H = He® le + Ic® He + Hint For analysis, we 
can choose Lindblad master equation^' 

p = ^[H,p] ~Y.^n{LlLnP + pL+Ln - 2LnpL+), (19) 

n 

which preserves the positivity of density matrix and it is Markovian but it is not 
general form of exact master equation. Other choice is Wigner transform of master 
equation: 

W = {H, W}pB + (20) 
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and it is more preferable for us. In the next Section we consider the variational- 
wavelet approach Refs. [2H7| for the solution of all these Wigner-like equations (2), 
(8), (13), (20) for the case of an arbitrary polynomial U{q,p), which corresponds 
to a finite number of terms in the series expansion in (13), (20) or to proper finite 
order of h. Analogous approach can be applied to classical counterpart (5) also. 
Roughly speaking, wavelet analysiJ^ is some set of mathematical methods, which 
gives the possibility to take into account high-localized states, control convergence 
of any type of expansions and gives maximum sparse forms for the general type 
of operators in such localized bases. These bases are the natural generalization 
of standard coherent, squeezed, thermal squeezed states which correspond to 
quadratic systems (pure linear dynamics) with Gaussian Wigner functions. The 
representations of underlying symmetry group (afline group in the simplest case) 
on the proper functional space of states generate the exact multiscale expansions 
which allow to control contributions to the final result from each scale of resolution 
from the whole underlying infinite scale of spaces. 

5. Variational Multiresolution Representation 

5.1. Multiscale Decomposition for Space of States: Functional 
Realization and Metric Structure 

We obtain our multiscale/multiresolution representations for solutions of Wigner- 
like equations via a variational- wavelet approach. We represent the solutions as 
decomposition into localized eigenmodes (regarding action of affine group, i.e. hid- 
den symmetry of the underlying functional space of states) related to the hidden 
underlying set of scales: 



where value ic corresponds to the coarsest level of resolution c or to the internal scale 
with the number c in the full multiresolution decomposition (MRA) of the underly- 
ing functional space (L^, e.g.) corresponding to the problem under consideration^ 



and p = (pi,p2,---)i Q = (9i,92,---): = (Pi, 9ii • ■ • .Pi, ft) are coordinates 
in phase space. In the following we may consider as fixed as variable numbers of 
particles. We introduce the Fock-like space structure (in addition to the standard 
one, if we consider second-quantized case) on the whole space of internal hidden 
scales 



oo 



W^it,q,p) = ^W^it,q,p), 



(21) 



Vc C Vc+l C Vc+2 C . . . 



(22) 



n 



(23) 



n 




(24) 
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where Wp{xi,. ..,Xp;t) e HP, H° = C, HP ^ L^{R^p) (or any different proper 
functional space), with the natural Pock space like norm: 

{W,W) = WS+Y^ I Wf{xl,...,x,■t)\{^l,. (25) 

i i=l 

First of all, we consider W — W{t) as a function of time only, W E L^{R), via 
multiresolution decomposition which naturally and efficiently introduces the infinite 
sequence of the underlying hidden scales We have the contribution to the final 
result from each scale of resolution from the whole infinite scale of spaces (22). The 
closed subspace Vj{j G Z) corresponds to the level j of resolution, or to the scale j 
and satisfies the following properties: let Dj be the orthonormal complement of Vj 
with respect to Vj+i: Vj+i = Vj ^Dj. Then we have the following decomposition: 

oo 

{Wm= D, = V^D„ (26) 

— C)0<j'<oo i— 

in case when Vc is the coarsest scale of resolution. The subgroup of translations 
generates a basis for the fixed scale number: spanj.g2{2^/^5'(2^i — k)} — Dj. The 
whole basis is generated by action of the full afhne group: 

spanfeez,,ez{2^/2^(2^i - k)} = span,^^.g^{vl/,. = {W{t)}. (27) 
5.2. Tensor Product Structure 

Let sequence {Vj},Vj C L^{R) correspond to multiresolution analysis on time axis 
and {Vj^'}, C L^{R) correspond to multiresolution analysis for coordinate 

Xi, then Vp^^ = Vj^ ® . . .® 1/^^" ® Vj corresponds to multiresolution analysis for 
n-particle distribution function Wn{xi, . . . , Xn] t). E.g., for n = 2: 

V^ = {/ : /(xi,X2) = Yl <^k,M^\xi - ki,X2 - fe), ak„k, G ^'(^')}, (28) 

ki ,^2 

where <j)'^{xi,X2) = (f>^{xi)(f>'^{x2) — (t)^ ® (f>'^{xi,X2), and (l)^{xi) = (f>{xi) form a 
multiresolution basis corresponding to {V^^'}. If {(j)^{xi — £)}, £ G Z form an or- 
thonormal set, then 4''^{xi — ki,X2 ~ ^2) form an orthonormal basis for Vq . Action 
of affine group provides us by multiresolution representation of L^{R^). After intro- 
ducing detail spaces D^, we have, e.g. V^ = V^ © D^. Then 3-component basis for 
Dq is generated by translations of three functions 

*2 ^ 0l(a;i) ® ^'2(^2), 'i'l = "f^Xi) ® 02(X2), = *i(xi) ® ^\X2). (29) 

Also, we may use the rectangle lattice of scales and one-dimensional wavelet decom- 
position : 

f{xi,X2)=Y^ «)%fc > %£(g)*j,fc(a;i,a;2), (30) 

ij;j,k 

where bases functions ^P^^f (g) 'ifj^k depend on two scales 2~* and 2"-' . 
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After construction the multidimensional bases we obtain our multiscale/multi- 
resolution representations for observables (symbols), states, partitions via the vari- 
ational approaches in Refs.[2Hl|as for c-BBGKY as for its quantum counterpart and 
related reductions but before we need to construct reasonable multiscale decompo- 
sition for all operators included in the set-up. 



5.3. FWT Decomposition for Observables 

One of the key point of wavelet analysis approach, the so called Fast Wavelet Trans- 
form (FWT) l^, demonstrates that for the large classes of operators the wavelet-like 
functions are best approximation for true eigenvectors and the corresponding ma- 
trices are almost diagonal. So, powerful FWT provides the maximum sparse form 
for different classes of operatorsl^l. Let us denote our (integral/differential) operator 
from equations under consideration as T (L^(i?") — > L'^(R")) and its kernel as K. 
We have the following representation: 

<T/,5>=y j K{x,v)f{v)g{x)dxAy. (31) 

In case when / and g are wavelets Lpj^k = 2^/'^ip{2^ x — k), (21) provides the standard 
representation for operator T. Let us consider multiresolution representation . . . C 
V2 C Vi C Vo C V-i C V-2 ■ ■ ■■ The basis in each Vj is (pj_k{x), where indices 
k,j represent translations and scaling respectively. Let Pj : L'^{R'^) Vj {j E Z) 
be projection operators on the subspace Vj corresponding to level j of resolution: 
{Pjf){x) = J2k ^ fj'Pj,k > ^3,k{x). Let Qj = Pj-i — Pj be the projection operator 
on the subspace Dj (V,-i — Vj ® Dj), then we have the following representation 
of operator T which takes into account contributions from each level of resolution 
from different scales starting with the coarsest and ending to the finest scales^ 

(QoTQj + QjTPj + P.TQj). (32) 

We need to remember that this is a result of presence of affine group inside this con- 
struction. The non-standard form of operator representation is a representation of 
operator T as a chain of triples T — {Aj^ Bj, Fjjjgz, acting on the subspaces Vj and 
Dj-. Aj : Dj Dj,Bj : Vj — > Dj,Tj : Dj — > Vj, where operators {Aj,Bj,Tj}j^z 
are defined as Aj — QjTQj, Bj — QjTPj, Tj = PjTQj. The operator T admits 
a recursive definition via 

_ / Aj+i B.j+1 



where Tj ~ PjTPj and Tj acts on Vj : Vj ^ Vj. So, it is possible to provide the 
following "sparse" action of operator Tj on elements / of functional realization of 
our space of states H: 

iTjf){x) = ^ [ 2-^' ^r,/,- ) fc(x), (34) 
kez \ e J 
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in the wavelet basis (pj,k{x) — 2 ^^^(p{2 — k), where 

f,,k-i = 2-^/2 J f{x)ip{2-^x -k + £)dx (35) 

are wavelet coefficients and are the roots of some additional linear system of 
equations related to the "type of localization" So, we have the simple linear 
parametrization of matrix representation of our operators in localized wavelet bases 
and of the action of this operator on arbitrary vector/state in proper functional 
space. 



5.4. Variational Approach 

Now, after preliminary work with (functional) spaces, states and operators, we may 
apply our variational approach from [2]-[7]. Let L be an arbitrary (non)linear diffe- 
rential/integral operator with matrix dimension d (finite or infinite), which acts on 

some set of functions from L^(il®"): 'J = 'I'(t, 2:1, X2, • • •) = {^^i^j xi, X2, ■■■),■■ ■, 
'i!'^(t, xi, X2, • • ■)) 1 e C R®, n is the number of particles: 

L'f = L{Q,t,Xi)-i'{t,Xi)^0, (36) 
Q = Qdo.di.d2,...{t,xi,X2, ■ ■ . ,d/dt,d/dxi,d/dx2, . . . , J Hk) 

= E ''^on^...Xt,x,,X2,. ..){-)"{— )\—)\. 



Let us consider now the N mode approximation for the solution as the following 
ansatz: 



N 



'b^{t,xi,X2, ■ ■ ■) = E aiaiii2...Ai„ Bi^ . . .{t,xi,X2, ■ ■ ■) (37) 



=1 



We will determine the expansion coefficients from the following conditions (related 
to proper choosing of variational approach): 



'koMM,...^ J {L^'')Ak„{t)Bkdxi)CkAx2)dtdxidx2... = 0. (38) 

Thus, we have exactly dA^" algebraical equations for dN" unknowns aig^i-^^.... This 
variational approach reduces the initial problem to the problem of solution of func- 
tional equations at the first stage and some algebraical problems at the second one. 
It allows to unify the multiresolution expansion with variational construction in 
Refs. 

As a result, the solution is parametrized by the solutions of two sets of reduced 
algebraical problems, one is linear or nonlinear (depending on the structure of the 
generic operator L) and the rest are linear problems related to the computation 
of the coefficients of reduced algebraic equations. It is also related to the choice of 
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exact measure of localization (including class of smoothness) which are proper for 
our set-up. These coefficients can be found by some functional/ algebraic methods 
by using the compactly supported wavelet basis functions or any other wavelet 
families'^. As a result the solution of the equations/hierarchies from Section 4, as 
in c- as in q-region, has the following multiscale or multiresolution decomposition 
via nonlinear high-localized eigenmodes 

W{t,xi,X2,...) = aijU' (g)V^t,xi,...), 

(ij)ez2 

V^t) - V^^^'^'^it) + J2 y/i^itl (39) 

1>N 

U^(xs) ~ Uj^j (xs) + ^ U^^{k'!^Xs), ~ 2™, 

m>Af 

which corresponds to the full multiresolution expansion in all underlying time/space 
scales. The formulae (39) give the expansion into a slow part and fast oscillating 
parts for arbitrary N, M . So, we may move from the coarse scales of resolution to 
the finest ones for obtaining more detailed information about the dynamical pro- 
cess. In this way one obtains contributions to the full solution from each scale of 
resolution or each time/space scale or from each nonlinear eigenmode. It should 
be noted that such representations give the best possible localization properties in 
the corresponding (phase)space/time coordinates. Formulae (39) do not use pertur- 
bation techniques or linearization procedures. Numerical calculations are based on 
compactly supported wavelets and wavelet packets and on evaluation of the accu- 
racy on the level N of the corresponding cut-off of the full system regarding norm 
(25): \\W^+^-W^\\ <e. 

6. Modeling of Patterns 

To summarize, the key points are: 

1. The ansatz-oriented choice of the (multidimensional) bases related to some 
polynomial tensor algebra. 

2. The choice of proper variational principle. A few projection/ Galerkin-like 
principles for constructing (weak) solutions can be considered. The advantages of 
formulations related to biorthogonal (wavelet) decomposition should be noted. 

3. The choice of bases functions in the scale spaces Dj from wavelet zoo. They 
correspond to high-localized (nonlinear) excitations, nontrivial local (stable) dis- 
tributions/fluctuations or "continuous qudits". Besides fast convergence properties 
it should be noted minimal complexity of all underlying calculations, especially in 
case of choice of wavelet packets which minimize Shannon entropy. 

4. Operator representations providing maximum sparse representations for ar- 
bitrary (pseudo) differential/ integral operators d//dx, d"//da;", / T{x,y)f{y)dy), 
etc. 
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5. (Multi)linearization. Besides the variation approach we can consider also a 
different method to deal with (polynomial) nonlinearities: para-products-like de- 
compositions. 

To classify the qualitative behaviour we apply standard methods from general 
control theory or really use the control. We will start from a priori unknown coefh- 
cients, the exact values of which will subsequently be recovered. Roughly speaking, 
we will fix only class of nonlinearity (polynomial in our case) which covers a broad 
variety of examples of possible truncation of the systems. As a simple model we 
choose band-triangular non-sparse matrices (a^ ). These matrices provide tensor 
structure of bases in (extended) phase space and are generated by the roots of the 
reduced variational (Galerkin-like) systems. As a second step we need to restore the 
coefficients from these matrices by which we may classify the types of behaviour. We 
start with the localized mode, which is a base mode/eigenfunction, which was con- 
structed as a tensor product of the two base functions. Fig. 8, 11 below demonstrate 
the result of summation of series (39) up to value of the dilation/scale parameter 
equal to four and six, respectively. It's done in the bases of symmlets [9] with the 
corresponding matrix elements equal to one. The size of matrix of "Fourier-wavelet 
coefficients" is 512x512. So, different possible distributions of the root values of the 
generic algebraical systems (38) provide qualitatively different types of behaviour. 
Generic algebraic system (38), Generalized Dispersion Relation (GDR), provide the 
possibility for algebraic control. The above choice provides us by a distribution 
with chaotic- like equidistribution. But, if we consider a band- like structure of ma- 
trix (fly ) with the band along the main diagonal with finite size (<C 512) and values, 
e.g. five, while the other values are equal to one, we obtain localization in a fixed 
finite area of the full phase space, i.e. almost all energy of the system is concen- 
trated in this small volume. This corresponds to waveleton states ^ and is shown 
in Fig. 9, constructed by means of Daubechies-based wavelet packets. Depending 
on the type of solution, such localization may be conserved during the whole time 
evolution (asymptotically-stable) or up to the needed value from the whole time 
scale (e.g. enough for plasma fusion/confinement in the case of fusion modeling by 
means of c-BBGKY hierarchy for dynamics of partitions). 

7. Conclusions 

By using wavelet bases with their best phase space localization properties, we can 
describe the localized (coherent) structures in quantum systems with complicated 
behaviour (Figs. 8, 11). The numerical simulation demonstrates the formation of 
different (stable) pattern or orbits generated by internal hidden symmetry from 
high-localized structures. Our (nonlinear) eigcnmodes are more realistic for the 
modeling of nonlinear classical/quantum dynamical process than the correspond- 
ing linear gaussian-like coherent states. Here we mention only the best convergence 
properties of the expansions based on wavelet packets, which realize the minimal 
Shannon entropy property and the exponential control of convergence of expansions 
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like (39) based on the norm (25). Fig. 9 corresponds to (possible) result of super- 
selection (einselection) [1] after decoherence process started from entangled state 
(Fig. 12); Fig. 10 and Fig. 13 demonstrate the steps of multiscalc resolution (or de- 
grees of interference) during modeling (quantum interaction/evolution) of entangled 
states leading to the growth of degree of entanglement. It should be noted that we 
can control the type of behaviour on the level of the reduced algebraical variational 
system, GDR (38). 

Let us finish with some phenomenological description which can be considered 
as an attempt of qualitative description of the quantum dynamics as a whole and 
in comparison with its classical counterpart. It is possible to take for reminis- 
cence the famous Dirac's phrase that "an electron can interact only itself via the 
process of quantum interference". Let G be the hidden/internal symmetry group 
on the spaces of quantum states which generates via MRA (22), (26) the multi- 
scale/multiresolution representation for all dynamical quantities, unified in object 
0{t), such as states, observables, partitions: 0^{t) = {t{j^{t),Op^{t),W^{t)}, where 
i is the proper scale index. Then, the following commutative diagram represents 
the details of quantum life from the point of view of representations of G on the 
chosen functional realization which leads to decomposition of the whole quantum 
evolution into the proper orbits or scales corresponding to the proper level of reso- 
lution. Morphisms W{t) describe Wigner-Weyl evolution in the algebra of symbols, 
while the processes of interactions with open World, such as the measurement or 
decoherence, correspond to morphisms (or even functors) m{t) which transform the 
infinite set of scales characterizing the quantum object into finite ones, sometimes 
consisting of one element (demolition/destructive measurement). 

W{t) 

{0\t,)} {0^it2)} 

i m{ti) i m{t2) 

Wit) 

{0^^{t,)} {0^=(t2)}, 

where reduced morphisms W{t) correspond to (semi) classical or quasiclassical evo- 
lution. So, qualitatively. 

Quantum Objects can be represented by an infinite or sufficiently large set of 
coexisting and interacting subsets like (22), (26) while 

Classical Objects can be described by one or few only levels of resolution with 
(almost) suppressed interscalc sclf-intcraction. It is possible to consider Wigner 
functions as some measure of the quantum character of the system: as soon as it 
becomes positive, we arrive to classical regime and so there is no need to consider 
the full hierarchy decomposition in the representation (21). 

So, Dirac's self-interference is nothing but the multiscale mixture/intermittency. 
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Certainly, the degree of this self-interaction leads to different qualitative types of be- 
haviour, such as localized quasiclassical states, separable, entangled, chaotic etc. At 
the same time the instantaneous quantum interaction or transmission of (quantum) 
information from Alice to Bob takes place not in the physical kinematical space- 
time but in Hilbert spaces of states in their proper functional realization where 
there is a different kinematic life. To describe a set of Quantum Objects we need 
to realize our Space of States (Hilbert space) not as one functional space but as 
the so-called and well known in mathematics scale of spaces, e.g. ^, ^ The 
proper multiscale decomposition for the scale of space provides us by the method 
of description of the set of quantum objects in case if the "size" of one Hilbert 
space of states is not enough to describe the complicated internal World. We will 
consider it elsewhere, while here we considered the one-scale case (to avoid possible 
misunderstanding we need to mention that one-scale case is also described by an 
infinite scale of spaces (26), but it is internal decomposition of the unique, attached 
to the problem, Hilbert space). 
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Fig. 2. Kick. 
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Fig. 5. MRA for Kick. 




Fig. 7. MRA for RW-fractal. 
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Fig. 8. Level 4 MRA. 




Fig. 9. Localized pattern, (waveleton) Wigner 
function. 
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Fig. 11. Level 6 MRA. 
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Fig. 12. Entangled- like Wigner function. 
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Fig. 10. Interference picture on the level 4 ap- Fig. 13. Interference picture on the level 6 ap- 
proximation for Wigner function. proximation for Wigner function. 



